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We examine the mathematical and physical significance of the spectral density (j{lu) introduced 
by Ford in defining the contribution of each frequency to the renormalised energy density of 
a quantum field. Firstly, by considering a simple example, we argue that a{ij) is well defined, 
in the sense of being regulator independent, despite an apparently regulator dependent definition. 
We then suggest that cr(a;) is a spectral distribution, rather than a function, which only produces 
physically meaningful results when integrated over a sufficiently large range of frequencies and with 
a high energy smooth enough regulator. Moreover, a-{ui) is seen to be simply the difference between 
the bare spectral density and the spectral density of the reference background. This interpretation 
yields a simple 'rule of thumb' to writing down a (formal) expression for o-{(jj) as shown in an 
explicit example. Finally, by considering an example in which the sign of the Casimir force varies, 
we show that the spectrum carries no manifest information about this sign; it can only be inferred 
by integrating o-(ctj). 

PACS numbers: ll.lO.Gh, 03.70.+k 



I. INTRODUCTION 

Can one assign a frequency spectrum to the finite 
Casimir energy? This question was first posed by Ford 
in [l[ where such frequency spectrum, henceforth dubbed 
Casimir spectrum, was considered for two examples: a 
massless scalar field in R x and R'^ x S^. The spec- 
trum is encoded in a spectral density <7{lu) which must 
obey two obvious criteria: i) integrated over all frequen- 
cies it must yield the renormalised vacuum energy 

1 f+°° 

Pren = (0| Too 1 0) ren = TT Cr{uj)duJ , (1.1) 

Jo 

where the quantum field has energy momentum tensor 
Tpi, in a space-time admitting a globally time-like Killing 
vector field d/d^o; ii) it should not depend on the renor- 
malisation/regularisation procedure. Whereas the first 
criterion is easy to check, the second one seems hard to 
prove. Indeed, the very definition of the spectral den- 
sity introduced in [l| was based on a specific rcgular- 
isation/renormalisation scheme, namely point splitting. 
This definition is 

/+CXJ 
(0|Too(r)|0).e„e''""dT , (1.2) 

where r is a point-splitting regulator along the time com- 
ponent. We shall show in Section IIIDl however, that 
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this definition, albeit apparently tied to point-splitting 
regularisation, is quite general, and thus, that in this re- 
spect the Casimir spectrum is well defined. From the 
same analysis we shall learn that renormalised spectral 
density 17(0;) is simply the difference between the bare 
spectral density and the spectral density of the refer- 
ence background. This provides us with an immediate 
tool to write down a formal expression for it. But if an 
appropriate regulator is not included, this expression is 
meaningless. Moreover, we shall see that even with an 
appropriate regulator, only the integration of this spec- 
tral density carries some physical information, and only 
if the integration extends through a sufficiently large fre- 
quency interval. Thus one should regard cr{uj) as a spec- 
tral distribution, rather than a function. 

As another test on the physical information carried by 
the Casimir spectrum, we analyse the case of a scalar 
field, with mass parameter and coupling to the curva- 
ture, in an Einstein Static Universe. It has been shown 
in Q, that the renormalised energy momentum ten- 
sor associated to the vacuum fiuctuations may obey or 
violate the strong energy condition, depending on the 
various parameters, therefore potentially sourcing repul- 
sive gravity. Our analysis shows, however, than no hint 
about the sign of the Casimir force may be obtained, 
in this case, from the Casimir spectrum. Thus, although 
we find that the Casimir spectrum is mathematically well 
defined, we find no evidence about it carrying more phys- 
ical information than that carried by the Casimir energy 
itself. 

This paper is organised as follows. In section 2 we 
start by reviewing the computation of the Casimir spec- 
trum for a massless scalar field in M x 5^, discussing the 
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generality of the definition (jl.2[) as well as its interpre- 
tation and physical content. In section 3 we discuss the 
Casimir spectrum on the Einstein Static Universe, where 
we show that no noticeable difference is seen when the 
gravitational effect of the Casimir energy varies from at- 
tractive to repulsive. In Appendix |X] we consider another 
example, a scalar field in R'^ x S'^, that confirms the ar- 
guments and interpretation provided in Section 2. 



II. RECONSIDERING R x 

We start by reconsidering the case of a massless scalar 
field on a 2-dimensional space-time with one periodic 
space-like direction, of length L. The space-time is then 
Rx S*^. We shall review the computation of the spectrum 
using a weight function [T], 0] and then compute it using 
point-splitting. 



A. Weight functions 

The regularised energy density (with a weight func- 
tion) is 

Pxo,m = ^ X! ^nW{uJn)xo,7n , (2.1) 
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with LUn = 2Ti\n\/L. We can apply the Abel-Plana for- 
mula as long as the condition 

lim e-'^''y\G{x±iy)\ =0 , 

i/->oo 

is satisfied (where G{n) is the function inside the sum). 
The choice of weight function in [i| is 



W{2ttx/L), 



\ 2m+l 



xo 



{2m) 
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where m £ N. This clearly obeys the conditions of appli- 
cability of the Abel-Plana formula. Applying the formula 
we get 
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ujW{uj)duj 
t {W{2TTit/L) + W{-2-Kit/L)) 
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dt . (2.3) 



The first term renormalises the cosmological constant 
since it is a background independent constant. The sec- 
ond term is identified as the spectral function. Note that 
the weight functions (|2.2p were chosen as to become delta 
functions 5{x — Xq) in the large m limit, which when inte- 
grated over xq give 1. Also, xq = Lujo/{2tt); the spectral 
function so defined is then normalised to give the correct 



energy density when integrated with this differential. Re- 
placing by the explicit form of the weight function we get 

, , 27r(-l)™+i /2m\'"+' 



dt . (2.4) 
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If we expand the denominator as a geometric series and 
perform the integrals, then the sums can be identified as 
Hurwitz zeta functions and the final expression coincides 
with that of [I: 



(-l)™+i(2m+ 1) / m \ 
[Xo) = 



2m + 1 



m \ m 
Cff 2to + 2, 1 - i + Cff 2m + 2, 1 + i 



(2.5) 



The plot in Fig. [T] reveals that in the limit of large 
m the spectrum becomes just the difference between lo- 
calised delta functions and a straight line with negative 
slope which matches the flat space spectrum. Indeed, 
re-writing the energy density as 



+00 



p=— a;„ = / ^^^^ ^i^-^n) 

n— — OQ n— — 00 

2 /•+00 +00 
= — / dL0Uj'S^5{uj — LOn) (2.6) 

and subtracting the flat space contribution (assuming a 
massless fleld), we get 
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+ 00 



dk 



2./-00 27r 



\ ■ (2.7) 



Thus, in this simple case, it is clear that the spectral 
density is really 

a{uj) = LO ^27r ^5{uj - a;„) - 1^ . (2.8) 

The terms in parenthesis show that the continuum fre- 
quencies have a negative uniform weight, whereas the 
discrete frequencies have a delta function uniform weight. 
The difference corresponds exactly to the limit of m — > 00 
in Fig. m 



B. Point Splitting 

The first form of the spectrum obtained in the previous 
section, depended on a particular choice of weight func- 
tion. But the same result is obtained with point splitting 
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FIG. 1. Left: Spectral function, am{xo) from weight functions for m = 1000, Right: Spectral function from weight functions 
for m = 10000. 



regularisation (or any other well defined regularisation). 
To show this, let us first define the point split energy 
density. As shown in the weight function is a cosine 
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iToo(r)lO) 
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a;„cos(a;„r) 



(2.9) 



Now, care must be take when applying the Abel-Plana 
formula. When the argument of the cosine is imagi- 
nary there is an exponentially growing part. Such term 
comes from the substitution e^*^'^"'^ ^-i2TT{x+iy)T 



which marginally violates the convergence condition for 
the Abel-Plana formula. We can cure this by adding a 
small imaginary part to r (which amounts to adding an 
exponentially decaying part to the cosine). If we multi- 
ply by a factor of the form e"'^'^", and expand the cosine 
in terms of exponentials, then we have 



+00 



(2.10) 



This gives a plus e prescription for one of the exponentials 
in the cosine and a minus e prescription for the other one. 
Applying now the Abel-Plana formula yields 
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(2.11) 



where we have discarded the term 



The definition of the spectral function (|1.2p yields, 
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dujuj cos{ujr)e 



(2.12) 



which amounts to performing the renormalisation of this 
quantity. For r > 1 the Abel-Plana formula does not 
converge, however we know that the original sum con- 
verges for all T, so if we obtain an analytic expression 
when assuming r < 1, it must hold for all t because of 
the analyticity of the original sum. If we expand the de- 
nominator using the geometric series for 1/(1 — e"^'^^) 
and integrate term by term, we get 
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(fc+ 1+11)2 (fc-l^)2 



(2.13) 
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(fc-l^)2 (fc-I±H)2 (fc+l_M)2 



dT . 
(2.14) 



The poles of the integrand are above or below the real 
axis for =Fze respectively. We can close the contour on the 
upper part of the complex plane, because of the exponen- 
tial which will give a vanishingly small contribution from 
the circle at infinity, so we reduce the problem to the cal- 
culation of residues. Then, after some manipulations we 
get 



+ 00 

cr(w) = 2w^e"'"cos(wA:L) 
fc=i 



(2.15) 
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Note that keeping the convergence factor with the e is 
crucial to induce the correct contour in the integration. 
We can integrate the spectral function p.l5|) to check 
that we recover the Casimir energy 



•'^ fc=i 
~ 7rL2 ^ (e2 + fc2) 



' cos{ujLk) 



k=l 



6X2 



(2.16) 



Using the Fourier series representation of the periodic 
repetition of the delta function one can show that the sum 
(|2.15p . with e 0, is equivalent to (|2.8p . Thus, using 
point-splitting we recover the spectral function obtained 
in the previous subsection using a weight function. 



C. Spectral function or spectral distribution? 

The natural question is: does o-(a;) carry some extra 
physical content, as compared to the energy density? Or 
in other words: is the spectral function a {to) meaningful 
only when integrated over the whole frequency range, or 
can one assign a physical value to an integration over an 
interval, J^"'^'^^ a{uj)duj ? 

To answer this question we can integrate cr(a;) in (|2.15p 
convoluted with some other functions. For example we 
can take functions which are non-zero only in regions 
with only one delta function contribution. For example 
if we choose, for j G N 



1 , (27rj - 7r)/X < a; < (27rj + tt) /L 
, otherwise ' 

(2.17) 

then, using (|2.15p . 



fj(uj)a{uj) — 
Zn 







(2.18) 



This is exactly what we obtain from integrating fj{uj) 
with the delta distribution centred at 2'nj/L, and sub- 
tracting the result from integrating fj{u}) with w/(27r), 
cf . (|2.8p . But integration over the first (remaining) inter- 
val yields 



TT 

4X2 



(2.19) 



rather than —tt/GL'^, which is the renormalised energy 
density. This is a consequence of manipulating divergent 
infinite sums without regularising them. Of course, if 
one keeps the convergence factor e^^" in (|2.16p . and in- 
tegrates again over the full range of frequencies, one will 
obtain the correct result. This makes clear that the spec- 
tral function a{u;) is only meaningful when regularised. 

What if we include a regulator in cr(a;)? Does the spec- 
tral function then acquire a physical meaning when inte- 
grated over a compact range of frequencies (as opposed 



to all frequencies)? A hint to answer this question comes 
from work on semi-transparent boundaries [sj. Therein, 
the Casimir force is given by a difference of pressures in- 
side the plates and outside the plates, computed from 
scattering of modes inside and outside the plates. In 
such computation there is a physical coefficient playing 
the role of a reflection factor, which becomes the physical 
regulator. A related idea was discussed in Q. 

Assuming that the spectral function has a physical 
meaning, we could think about integrating it over a fi- 
nite range of frequencies, to mimic the situation of trans- 
parency in the remaining range of frequencies. Let us first 
consider this is done in a sharp way, i.e. by introducing 
a Hcaviside step function H{x) cut-off: 



Fi(g,wo) = J —(7{uj)H{q- 



(2.20) 



In Fig. [2] we plot Fi(q,ajo). It oscillates wildly and it 
does not converge to the renormalised energy density in 
the limit wo — >■ oo. 

However, we can introduce a smoother regulator, as 
we would expect in an actual physical situation such as 
in @. Using the exponential cut-off e^*^", so that we can 
interpret g = 1/e as the scale at which the theory is UV 
completed; then consider 



F2(q,c.o)^^""^a(c.)e-/« 



(2.21) 



This function is plotted in Fig. |3l It can be seen that 
the integrated value of the spectral function now ap- 
proaches asymptotically the Casimir energy density of 
the system. Thus, for the spectral function to be phys- 
ical, a regulator that suppresses the UV frequencies is 
crucial and it should correspond to an actual physical 
cut-off of the theory. It is a remarkable fact that the 
Casimir force does not depend on the details of such regu- 
lator, but only on the existence of a high energy "smooth 
enough" damping of the modes. 



D. Interpretation of (T(a;) and computational 
shortcut 

Inspired by the previous example, and since the defini- 
tion of the spectral density (|1.2p is not derived from first 
principles, let us reconsider it starting from the (bare) 
energy density formula. Keeping the convergence fac- 
tor e-"'^'^ 
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/ '^tY. (^('^ -^n) + 5{u; + iJn)) e-^"" .(2.22) 
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FIG. 2. Integral of the spectral function regularised with a step function at q, up to an energy uq Left: _Fi (100, cjq), Right: 
i^i(1000,a;o). 
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FIG. 3. Integral up to an energy cjq, of the spectral function regularised with an exponentially decaying function Left: 
-f2(100, cjo), Right: ^2(1000, tJo). The asymptotic result is indicated by the horizontal line in the zoomed plots on the right. 



Observe that the second delta function does not con- 
tribute. From this we see that the bare spectral function 
with the normalisation is 



27r 



+ 00 



fba 



{5{UJ - LOn) + 5{uJ + U)ri)) & 



n=l 



(2.23) 



If we replace here the Fourier expansion of the delta func- 
tion 



-t-oo 
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(0|Too(r)|0)e'™dr . 



(2.24) 
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So the fact that the spectral function is the Fourier trans- 
form of the point spht energy density is not restricted to 
using point splitting regularisation; indeed we have not 
used it in this argument. This expression appears simply 
from Fourier expanding the delta functions. 

The lesson from this example is that the renormalised 
spectral density cr(w) is simply the difference between the 
bare spectral density above and the spectral density of 
the reference background. In Appendix A we show that 
this simple 'rule of thumb' yields the correct result for 
the computation of the spectral density in another well 
known example: a massless scalar field in M'^ x S'"'^. 



III. APPLICATION TO R X 



the global timelike Killing vector). This can be done by 
solving the classical equations of motion for the scalar 
field, then quantising the theory and finally computing 
expectation values of observables in such a state. Since 
we are interested in the energy density we have to look 
at the expectation value of the energy momentum tensor 



2e - 2 ) 



(3.3) 



This procedure was explained in detail in 0, [1] , so here 
we simply summarise the main results needed for our 
study in 3 + 1 dimensions. The scalar field is expanded 
in the Heisenberg picture using the eigenmodes 



We will now ask another question concerning the phys- 
ical content of the spectral function: does it carry any 
hint concerning the sign of the Casimir force? It is well 
known that the Casimir force can either be attractive or 
repulsive depending on details of the geometry, topology 
or boundary conditions that originates it (see [7[ for a 
review). What determines the sign of the Casimir force 
is still an open problem. 

This question is particularly relevant since Casimir 
forces are a dominant interaction at the nanometer scale 
and therefore of relevance for the construction of nano- 
devices. Repulsive Casimir forces (see e.g. [1,0) could 
allow quantum levitation and therefore avoid collapse 
and permanent adhesion between nearby surfaces, which 
is a problem for microelectromechanical systems. More- 
over, repulsive Casimir forces are not just a theoretical 
possibility since there is some experimental evidence [1- 

In order to tackle this question we shall consider an 
example where it has been shown that the effect of the 
Casimir energy can change from attractive to repulsive, 
when one varies the parameters of the system: a scalar 
field with a mass parameter and a coupling to the Ricci 
scalar in the Einstein Static Universe (ESU) d, Q . 



A. Setup 

We consider a scalar field theory in the 3 -t- 1 dimen- 
sional ESU 



ds%gjj = -dt^ + R^dUss. . 



(3.1) 



The Ricci scalar of the geometry is 7?. = 6/i?^, where the 
radius of the universe is R. This geometry is perturbed 
by the quantum vacuum fluctuations of a real scalar field 
with action 



(3.2) 

We are interested in analysing the energy spectrum of 
the quantum field in a vacuum state (with respect to 



(t)a{t,x) 



(3.4) 



where a = {i, mn.mL} : ^ G No A — £ < mL,rnii < £ £ Z, 
Ya{x) are the hyperspherical harmonics and x denotes 
the angular coordinates on S^. The eigenfrequencies are 



with 



, ^+1 2_ 2 

k = — ^ , a = fi 



6^-1 



(3.5) 



(3.6) 



The Vacuum Expectation Value (VEV) of any compo- 
nent of T*^ can be expressed in terms of (0| (T*)°g |0) 
using the conservation of energy. Its bare value is for- 
mally divergent, so it requires a regularisation procedure 
followed by renormalisation which we address in the next 
section. 



B. The spectrum on the ESU from a generic 
regularisation 

As seen in the renormalisation of the energy mo- 
mentum tensor on the ESU can be performed just by 
specifying some generic physical properties for the UV 
completion of the theory, encoded in a generic regulator. 
Then the regularised energy density is simply 



^ +00 



(3.7) 
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where V is the volume of 5''^, g is a generic regularisation 
function that goes to 1 in the infrared, 7 ^ 1 is the 
regularisation parameter, L is the typical length scale of 
the system, and the degeneracy of the modes is 



(3.8) 



Using the Abel-Plana formula po can be re-written as 

P0 = -^ dt^yuP^^^Uj'^gi'^Luj) + Pren , (3.9) 
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where LOmm = a if > 0, otherwise LOmin = 0. The first 
integral is the quantity that is discarded and contains 
all the divergences that renormalise various gravitational 
couplings. In fact, to obtain the spectrum, this is exactly 
the quantity that we need to subtract from the discrete 
spectrum. This can be written as 



(3.10) 

So now, using our 'rule of thumb' we can just subtract 
the integrand from the discrete spectrum to obtain the 
spectrum 



(^ESU 



V 



+ 00 



1=0 



—R^uj\/ LiS^ — a^9{uj — oJnmi) gilLuj) . (3.11) 

Once again we can define the integral of the spectrum 
up to an energy loq and make the plot for various a. If 
for example we take the massless case, then we can show 
that for this case, the quantity p + 3p, which determines 
if the strong energy condition is obeyed, simplifies to 



p + 3p=2p , 



(3.12) 



so that it is simply the sign of p that determines if the 
gravitational effect of the Casimir energy is attractive 
(p > 0) or repulsive {p < 0). The change from a repulsive 
to an attractive effect can be obtained by varying a^R'^, 
or cquivalcntly ^ @j 01 ■ convenience wc define the 
integrated spectrum 



(3.13) 



with g{'^Luj) 



Some examples are exhibited in 



Fig. m for which one varies such that the effect varies 
from attractive to repulsive. The conclusion is that there 
is no manifest qualitative change in the spectrum and 
the variation in character of the Casimir force can only 
be observed upon integration of the spectral function. 
Thus, it seems that again, the Casimir spectrum carries 
no added value, relatively to the Casimir energy density, 
in terms of physical content. 



IV. CONCLUDING REMARKS 

In this paper we have reconsidered the Casimir spec- 
trum, defined in [ij. Based on the study of three dif- 
ferent examples (two of which had already been consid- 
ered in [l|) we have reached the following conclusions: 

1) despite the apparent regularisation scheme dependent 
definition (|1.2p . the Casimir spectrum is mathematically 
well defined, in the sense of being scheme independent; 

2) a simple interpretation for the spectrum is provided: 



the renormalised spectral density cr(u;) is simply the dif- 
ference between the bare spectral density and the spec- 
tral density of the reference background. This provides 
us with a simple rule to write down a formal expression 
for the spectrum; 3) the spectrum should, however, be 
faced as a distribution, rather than a function. That is, 
it is not differentiablc but it is intcgrable, and also it is 
only physically meaningful when integrated over a suf- 
ficiently large range of frequencies, with a high energy 
smooth enough regulator. Thus we find no evidence that 
the spectrum carries an added value, in terms of phys- 
ical content, relatively to the energy density itself. In 
particular, by considering an example where the effect 
of the Casimir energy can vary its character from repul- 
sive to attractive, we found that there is no manifest 
imprint of this change of character in the properties of 
the spectrum. We note that a possible relation between 
the Casimir spectrum and the sign of the Casimir force 
was previously discussed in 
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APPENDICES 
Appendix A: Check for 5*^ x 

In [l| the spectrum of the Casimir effect, for a massless 
scalar field in fiat spacetime, with one periodic spatial 
direction was computed using Green's functions. In this 
appendix we show that the shortcut suggested in Section 
Hid] yields the same result, which is also checked by using 
point-splitting regularisation. 

In this case, due to the boundary conditions, the field 
expansion is 

I (fkY^\al{k)(|)i{k,x^')+a\{k)(j3l{k,x^ 
•I f 

(Al) 

where the creation and annihilation operators obey stan- 
dard harmonic oscillator commutation relations, and 
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FIG. 4. Integral up to energy loq, of the spectral function regularised with an exponentially decaying function, for various 
choices of a^. The asymptotic result is indicated by the horizontal line in the zoomed plots (right). Note that — is the value 
for which the renormalised energy density vanishes; which can be seen on the corresponding zoomed plot up to an error of 10~*. 
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The first step is to eompute ttie vacuum expectation value of the regularised energy density from the energy 

momentum tensor 



|roo(.T-x')|o) = i(o|ao$ao'$' + iao'$'a"$ + xox'|o) 
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+ 00 
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J{ky{y-y')+k^{z-z')} 
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where we have used the field expansion and the canonical 
commutation relations for the creation and annihilation 
operators. 

Taking time separation only r = t — t', and introducing 
a convergence factor to make sure the sum converges we 
have 



1 



|roo(T)|o) = — 



2L J (27r) 



■ E 

l=-oo 



ujpk cos (Wf.fcTj e 



(A5) 

Before applying the Abel-Plana formula let us use the 
insight of Scction lll DI to anticipate the answer (since the 
spectrum was seen as a difference between the bare spec- 
tral density and the flat space spectral density) . The bare 
spectral density is obtained by inserting a delta function 
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So 



1 f d'l 



d'k 



^ e=-oD 

^ +00 ^ + oo 

^ E / dkkuJi^kH^} - uji,k) 







1^ ( l + 2^0(c.L-27r| 



2 / +°° 



The flat limit spectrum is clearly 



(A8) 



so we expect the spectrum to be the difference be- 
tween (1X71) and 



2ttL 



TT — ujL + 2n 



+ 00 

E' 



{ujL-27r\e\) . (A9) 



The function in parenthesis is just a (periodic) saw func- 
tion, which is a repetition of the first linear segment 
TT — u)L. This is exactly the same as the result in [l[ (ex- 
cept that the result therein has an incorrect overall minus 
sign). Note that wc have not included a converging fac- 
tor e^*^", which should be implicit for this distribution to 
make sense. 

Once again this result can also be obtained by applying 
the Abel-Plana formula to (jASP and performing manip- 
ulations similar to those in the appendix B of [S!]. Start 
with the point split energy density 



1 f d'^k 

P- = ^J E Wf.fcCOs(^,,fer)e-^'^^.'= , (AlO) 



(A7) and apply the Abel-Plana formula to get 



p,s(0|Too(r)|0)= J ^fccos(fcT)e-^^+^ J d\ J^"° dx J^ - cosh (27rV^^^^) cos (2^^^^^^) 

(All) 

I 



The first (divergent) term is a background independent 
constant, so it renormalises the cosmological constant. 
Since we are in flat space and all other curvature invari- 
ants are zero it can be shown that no further couplings 



get renormalised. The renormalised point split energy 
density is then (we have performed the angular integra- 
tion in d^q) 



10 



27r2 



+ 00 



dqq 



(A12) 



If wc expand the denominator in a geometric series once which we can check, by considering the Fourier series 
again, and exchange the order of integration, then we can of a saw-tooth function, to be exactly the same as the 



perform the integrals explicitly to give 



distribution (|A9p anticipated above - Fig. [5] 
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(A13) 



Once again if we take the Fourier transform, we can 
perform the integral by closing in the upper complex 
semi-plane and reduce the problem to the calculation of 
residues. The result is then 
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27ri 
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(A14) 




FIG. 5. Spectral function from point splitting for x R'^ 
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